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Abstract 



The iJL-T symmetry imposed on the neutrino mass matrix in the flavour basis is known to be 
quite predictive. We integrate this very specific neutrino symmetry into a more general frame- 
work based on the supersymmetric S'O(IO) grand unified theory. As in several other models, the 
fermion mass spectrum is determined by Hermitian mass matrices resulting from the renormaliz- 
able Yukawa couplings of the 16-plet of fermions with the Higgs fields transforming as 10, 126, 120 
representations of the 50(10) group. The ^-t symmetry is spontaneously broken through the 
120-plet. Consequences of this scheme are considered for fermion masses using both type-I and 
type-II seesaw mechanism. This scenario is shown to lead to a generalized CP invariance of the 
mass matrices and vanishing CP violating phases if the Yukawa couplings are invariant under the 
jjL-T symmetry. Small explicit breaking of the /i-r symmetry is then shown to provide a very good 
understanding of all the fermion masses and mixing. Detailed fits to the fermion spectrum are 
presented in several scenarios. One obtains a very good fit to all observables in the context of the 
type-I seesaw mechanism but type-II seesaw model also provides a good description except for the 
overall scale of the neutrino masses. Three major predictions on the leptonic mixing parameters in 
the type-I seesaw case are (1) the atmospheric mixing angle close to maximal, (2) close to 
the present upper bound and (3) negative but very small Dirac CP violating phase in the neutrino 
oscillations. 
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I. INTRODUCTION 



There exist variety of theoretical frameworks/specific models [l| which try to account for 
the large atmospheric mixing angle observed more than a decade ago. One class of theories 
attribute the maximal atmospheric mixing to the presence of some underlying flavour sym- 
metry. This would be a preferred alternative if the deviation of the atmospheric mixing angle 
from maximality is constrained to be very small. The simplest of such flavour symmetries is 
the fi-T symmetry which exchanges the mu and tau fields. This symmetry comes 

with an additional prediction that one of the three leptonic mixing angles namely, ^^^3 0] 
must be zero. 

fi-T symmetry is predictive and simple but it appears to have two shortcomings. Suc- 
cessful predictions follow only if it is an effective symmetry of the neutrino mass matrix 
in specific basis corresponding to a diagonal charged lepton mass matrix. The underlying 
flavour symmetry in general may not pick up this basis. Secondly, /i-r symmetry has been 
proposed with a view of explaining the mixing angles in the leptonic sector alone. It would 
be more desirable to have a symmetry providing overall understanding of complete fermionic 
mass spectrum. This can be done using the grand unified theory as the underlying frame- 
work. Various alternatives within such theories to simultaneously obtain small mixing in 
the quark sector and large mixing among leptons have already been proposed 

The renormalizable theories based on the 5*0(10) group are quite powerful in constraining 
the fermionic mass structures. The standard fermions are assigned to the 16 dimensional 
representation of the 5*0(10) group and they can obtain masses through symmetric couplings 
with 10 and 126 and antisymmetric couplings with the 120 dimensional representation of 
the Higgs fields. The minimal 5*0(10) model containing 10, 126, 126 and 210 representations 
has been extensively studied 0, 0, [lO, 11, 12, 0, Q, [l5|. In this model, the largeness of the 
atmospheric mixing angle gets related to the b — t Yukawa unification if neutrinos obtain 
their masses through the type II 16 1 seesaw mechanism [l3]. This interesting observation 
in [3| led to many detailed investigations 12, [isl, Q, [l5| which revealed the inadequacy of 
this simple picture. The supersymmetric version of the minimal model with type-II seesaw 
mechanism is constrained by two conflicting requirements. The overall neutrino mass scale 
is correctly reproduced in the model if the seesaw scale is about two to three orders of 
magnitude below the GUT scale. But the spectrum of the model in this case does not 
allow gauge coupling unification. Moreover, the type-II contribution to neutrino masses 
does not always dominate over the type-I contribution in the minimal model as would be 
required for the mechanism in [3] to go through. The conflict with the proton decay appears 
in the minimal model even if the neutrinos obtain their masses through the type-I seesaw 
These problems have led to studies of the non-minimal models containing an 



additional 120-plet of Higgs [l8|, [l9|, l20|]. Theoretical understanding of the largeness of the 



atmospheric mixing angle gets lost in all these approaches although one can choose the 
parameters to obtain the observed value. 

It would be welcome to integrate /x-r symmetry into the grand unified framework. This 
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has been done for the SU{5) model in 2l(]. We do so here in the more predictive 5*0(10) 
framework. There are several motivations for unifying 150(10) and /i-r symmetry. Rather 
than remaining a leptonic symmetry, such symmetry would provide constraining picture of 
both quark and lepton spectrum. Role of this symmetry in the description of the quark 
mixing is already discussed in 0, [sf. In addition, it can provide additional constraints and 
reduces the number of the Yukawa couplings which describe fermion masses. Some examples 



of models unifying 5*0(10) with other discrete symmetries can be found in 22 



We investigate the consequences of imposing a generalized fi-r symmetry exchanging the 
second and the third generation fields on a renormalizable 5*0(10) model. We deviate from 
the minimal model and add a 120 plet. This plays a crucial role in generating CP violation 
and the /i-r symmetry breakirig. Fermion masses with 120 plet have been discussed in several 
earher works H, 0, H H Q H H. Following H, H, 0, Q we impose also the 



, 119|, 120, 123|, 12J, 125|, 126|, 127J. Following [19|, C 
parity symmetry which leads to Hermitian mass matrices for all fermions thereby reducing 
the number of parameters compared to more general models. All the fermion masses and 
mixing are described in our approach in terms of 14 (15) real parameters in case of type-II 
(type-I) seesaw mechanism. They provide an excellent description of fermion masses and 
mixing in contrast to a general model employing 10 + 120 + 126 Higgs fields which needs 26 
31 parameters in the fermionic sector. Moreover, the (near) maximality of the atmospheric 
mixing and smallness of the angle ^'^^3 get related here to the approximately broken //-r 
symmetry. 

We define our model implementing fi-r symmetry and discuss its consequences in the next 
section. Section(III) presents numerical fits both in case of the type-II and type-I seesaw 
mechanism and discusses various predictions. Last section contains a summary. 



II. n-T SYMMETRIC 5*0(10) 

If /i-r symmetry is to be integrated with grand unification then a more general symmetry 
which exchanges the second and third generations of fermions should be imposed. Con- 
sequences of this generalization were first considered in 0]. It was subsequently noted [sl 
that this generalization automatically leads to understanding of why Cabibbo angle is larger 
than other two angles and a mild breaking of this symmetry was shown to lead to a correct 
description of the quark mixing angles and masses. Most of these works did not use the 
grand unified framework. Here, we consider a model based on the 5*0(10) (8> Z!^^^ eg) Z^^. The 



first Z2 corresponds to the generalized /i-r symmetry. The second Z2 symmetry called [26 1 
"parity" interchanges two components of the 16 field transforming as (4,2,1) and (4,1,2) 
under the Pati-Salam group decomposition of 5*0(10). 



Our basic formalism is similar to 19|, |20|, l26|, |27|. 16 dimensional fermions obtain their 



masses from coupling to three Higgs multiplets transforming as 10, 126 and 120 representa- 
tions under 5*0(10). The 5*0(10) breaking can be achieved with a 210-plet. An additional 
126-plet of Higgs is needed in the supersymmetric context to preserve the supersymmetry 
at the GUT breaking scale. These Higgs multiplets contain altogether six doublets with 
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quantum numbers of the minimal supersymmetric standard model (MSSM) field Hd and six 
with that of It is assumed that only two appropriate linear combinations of these Higgs 
doublets remain light and play the role of the and fields. This is achieved by the fine 
tuning conditions 1^. After this fine tuning, the resulting fermion masses can be written 

as 



- Cynass = flMffR + VlMd^R + ]^PlMlV'l + ^1^%MrUr + h.C. 



(1) 



where f = u,d,l denote the up and down quarks and the charged leptons respectively. The 



mass matrices appearing in the above equation can be suitably written (see 
details) as 

Md = H + F + iG , 
Mu = rH + sF + i t G , 
Ml = H -3F + i p a , 
Md = rH -3sF + i q G , 
Ml = tlF , 
Mr = r-^^F. 
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for 



(2) 



Here M£, denotes the neutrino Dirac mass matrix. Ml{Mr) is the Majorana mass matrix for 
the left (right) handed neutrinos which receives contribution only from the vacuum expecta- 
tion value (vev) of the 126 field. Gauge coupling unification in the minimal model requires 
that the vev contributing to Mr be close to the GUT scale. The dimensionless parameters 
r, s, t,p, q, tl and are determined by the CG coefficients, ratios of vevs and mixing among 
the Higgs fields 
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The matrices H, F, G originate from the fermion couplings to the 10, 126 and 120 fields 
respectively. (G) H, F are complex (anti) symmetric matrices in general. However, gener- 
alized parity makes them real. In addition, if all vevs and (hence r, s,t,p,q,rL,rR) are real 
then all the Dirac masses in eq.(l2]) are Hermitian and M^, Mr are real. 

We assume that the Higgs field in the 10 and 126 representations are invariant under 
the generalized /i-r symmetry while the 120 dimensional representation changes sign. This 
assumption allows spontaneous breaking of the fi-r symmetry. The resulting structures for 
H, F, G are given by 



H 



^ hii h 
\ 



12 

hi2 h22 h23 

hi2 h23 h 



\ 



22 y 



^ fu fl2 fl2 ^ 

/l2 /22 /23 
\ /l2 /23 /22 J 



G 



' gi2 -gi2 ^ 

-gi2 g23 
\ 9l2 -923 J 



(3) 



All the coefficients in these matrices are real. They satisfy 



S^{H,F,G)S = iH,F,~G) 



(4) 
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where 

/l o\ 

5=001 (5) 

\o 1 o) 

exchanges the second and the third generations. The effective neutrino mass matrix Aiu for 
the three hght neutrinos follows from eq. ([T]) and eq. : 

= tlF - trMdF-^MI = M'J + Mi . (6) 



Here tl^r are inversely related to the vev of the RH triplet component in 126. This vev may 
be identified with the GUT scale in the absence of any intermediate scale. In addition, they 
depend upon the details of the superpotential. Specific expressions for t^^r in the minimal 
case can be found in (isl . . The first term corresponds to the type-II seesaw while the 
second is the conventional type-I seesaw. In general, both contributions are present but one 
may dominate over the other. We shall be considering two separate cases corresponding to 
the type-II and type-I dominance respectively. 

The relations = f and = are major predictions and motivation for imposing the 
/i-r symmetry. These can arise if the effective neutrino mass matrix M.vf in flavour basis 
possesses a /i-r symmetry. Let us see how this can come about in our approach. It is easy 
to see that the fermionic mass matrices in our model satisfy: 

s-^MfS = m; , (7) 

S-'M'JS = M'J , (8) 
s-'Mis = Mi* . (9) 

f = u,d,l,D label the (Dirac) fermionic mass matrices. The Mi'^^ correspond to the type-I 
and II contributions to the light neutrino mass matrix, eq.Q. Let us note that 

• Eq.([S]) implies an exact /i-r symmetry for Mi^ . 

Ea s.([71P) correspond to an invariance under the generalized CP transformation defined 



a 
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as 

fa ^ ^SaP^'CJ^^ (10) 



If eq. (E]) represents the neutrino masses in flavour basis then one obtains the predictions 
= f and ^^3 = 0. 

If eq. ([9]) holds in the flavour basis then only the ^23 is maximal with definite correlations 
of with the CP violating phase Spmns 
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Ml is not diagonal here and hence these predictions do not follow immediately. It is 
still possible to recover these predictions even with a non-diagonal Mi. 
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Define 

UlMiUi = Di , (11) 

where Di is the diagonal mass matrix for the charged leptons. By factoring out a diagonal 
phase matrix Pi, the Ui can be written as: 

Ui = iliPi (12) 

The neutrino mass matrix in the flavour basis is then given by 

M.f = PjulMM^P; = PlM^fP; (13) 

The predictions of the fi-r symmetry are recovered if A4uf is /U-t invariant. This does not 
require a diagonal Mi. A general /i-r symmetric Ui satisfying S~^UiS = Ui will do the job 
in case of the type-II dominance. In case of the type-I dominance, one obtains 

S-'M,fS = Mlf 

provided Ui also satisfies the same equation. This makes it possible to recover the predictions 
of the fi-T symmetry for a non-diagonal Mi and obtain reasonably good fits to other fermion 
masses and mixin g. 

It is known [3I.I28I that with appropriate choice of Pi, Ui can be cast into the following 
form if Ml satisfies eq.([7]): 

^ Uii U21 U31 ^ 

Ui = wii W21 wsi , (14) 

V <i J 

with real uu. A unitary matrix with this form can be parametrized in terms of two angles 
and a phase. 



Cl S1C2 S1S2 

Ui = Pv\ ^ -^{ciC2-ies2) -^{ciS2 + iec2) 



(15) 



^ -73(ciC2 + ^es2) -j^{ciS2-iec2) J 



where e = ±1, Si^2 = sin 61^2 , Ci,2 = cos 61^2- C2 and S2 can be chosen positive with appropriate 
choice of Pi in eq. (fT2l) . 



iri\ 



Pr, = diag.(l, e e 

is a diagonal phase matrix. The above [// becomes /i-r symmetric if S2 = C2 and r] = 0. This 
defines a one parameter family of the leptonic mass matrices which lead to the prediction 
of the /i-r symmetry in case of the type-II dominance. We will use this form subsequently 
in our numerical analysis. 

There is an important but unwelcome feature associated with the generalized CP invari- 
ance of the mass matrices in eq.©. The CKM matrix in this case turns out to be real. To 
see this explicitly, we note that just as in case of Ui, the matrices Uu^d diagonalizing the up 
and down quark masses can be written as Uu^Pu^- Uu,d have the same form as the RHS of 
eq. (fT4|) with the replacement of uu with Uiu^id and wu with Wiu^id- The phase matrices Pu^ 
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can be absorbed in redefining the quark fields and the remaining part of the CKM matrix 
is given by 

Vij = {UlUdjij = UiuUjd + 2Re{wiuW*d) 

which is real since Uiu^id are real. 

One can generate CP violation in the model by breaking the generalized CP invariance of 
the mass matrices. This can be done in two ways. Either one allows complex vev for some 
of the Higgs doublets as in [l9| or one retains the real vev but allows breaking of the fi-T 
symmetry in the Yukawa couplings. In the following, we will discuss the second alternative. 



III. FITTING FERMION SPECTRUM WITH AND WITHOUT THE /x-r SYM- 
METRY 

We now discuss the numerical implications of our model in detail. We assume that either 
the type-I or the type-II term in the neutrino mass matrix dominates and carry out analysis 
separately in each of these two cases. Our input parameters are r, s, t,p, q, eq.([2]) , the real 
elements of the matrices G,H,F, eq.Q and the overall scales r^j L, eq.([6]). Parameter q is 
absent in the type-II case. An overall rotation R on G,H,F: {G,H,F) — > R^{G,H,F)R 
amounts to a choice of initial basis for the 16-plet of fermions. We can use this freedom to 
set say, = 0. This is done with a specific choice R = R23{j) Rui^u) R23{j) ■ Here Rij{0) 
denotes rotation in the ij^^ plane by an angle 6 and 



tan 26'^. 



2V2hu 



^"^ hii — h22 — h23 

This rotation amounts to redefinition of elements of F and G which still retain the same form 
as in eq.Q. We continue to use the same notation for the parameters of the redefined F, G. 
With the choice hi2 = 0, we have 14 (15) input parameters in case of type-II (type-I) seesaw 
dominance. These input parameters together generate 12 fermion masses and six mixing 
angles. As already remarked, the exact /i-r symmetric H, F, G are not able to generate CP 
violation. We introduce this CP violation by adding a small fi-r breaking difference between 
the 22 and 33 elements in H. This one additional parameter now leads to four CP violating 
phases, one in the CKM matrix and three in the PMNS matrix. 



Our choice of the values of the physical observables is based on numbers given in 13|, |20 
We reproduce them here in Table @ for convenience. 

The given numbers for quark masses and mixing correspond to the respective values at 
the GUT scale obtained from low energy values using MSSM and tan (3 = 10. The neutrino 
masses and mixing that we use are the low scale values but the effects of the evolution to 
Mgut on the ratio of the solar to atmospheric mass scale and on the mixing angles are known 
to be small for the normal hierarchical spectrum that we obtain here . While fitting, we 
omit the parameters r^, which define the overall scales of neutrino masses in case of the 
type-I and type-II seesaw respectively. The ratio of the solar and atmospheric mass scales 
and neutrino mixing parameters are independent of these overall scales and are used in our 
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"I'd 


1.03 ±0.41 


SOL 


(7.9 ±0.3) X 10-5 


rus 


19.6 ±5.2 




(2.2+Oi?) X 10-3 


rrih 


1063.6+^f,^ 

—86.5 


sin 


0.2243 ± 0.0016 


■niu 


0.45 ±0.15 


sin ^23 


0.0351 ±0.0013 


rric 


210.3273+i?-^^^^ 


sin 0?q 


0.0032 ± 0.0005 


mt 


on/foQ Q+30267.6 


sin2 0^2 


0.31 ± 0.025 


rUe 


0.3585 ± 0.0003 


sin2 


0.5 ±0.065 




75 671 '^+0-0578 


sin^ 


< 0.0155 


nir 


1292. 2l}| 


ScKM 


60° ± 14° 



TABLE I: Input values for quark and leptonic masses and mixing ang les at Mgut = 2 x 10^*' GeV 
and tan/3 = 10 which we use in our numerical analysis. 

definition of function instead of the individual neutrino masses. In addition, we assume 
Am^^^ to be positive corresponding to the normal neutrino mass hierarchy. Parameters 
tr^tl are fixed subsequent to minimization using the atmospheric scale. 

A. Numerical analysis: type-II seesaw 

We perform the minimization in three physically different cases. 

(A) In this case, we impose the conditions ^23 = f ^ind 6*5^3 = using a //-r symmetric 
Ui. As discussed in the earlier section, this is done using parametrization in eq. (fT5|) with 
S2 = C2 = The charged lepton mass matrix is then determined completely in terms 
of three masses and the angle 9i. Using the third of eq.([2]), the real and imaginary parts 
of Ml can be used to determine respectively elements of H in terms of that of F and 
elements of G in terms of p, the charged lepton masses and 61. fu also gets determined 
in terms of these parameters because of the choice hi2 = 0. Thus f22, f23, f 11, r, s,t,p, 61 
are the only free parameters which determine the 11 remaining observables- six quark 
masses, three angles of the CKM matrix, the solar angle and the solar to atmospheric 
mass ratio. The we minimize is defined in terms of these observables using the values 
and errors given in Table ([T]). The result of the minimization are shown in Tables fllTlIIip . 
One obtains a reasonably good fit to all observables except the down and bottom quark 
masses which are respectively ~ 1.5 and ~ 2.5 sigma away from their respective mean val- 
ues. All other observables are reproduced correctly with very small pulls as seen in the table. 

(B) In this case, we do not impose the maximality of 6*23 but include sin^ ^23 the to 
be minimized, sin^ 6*5^3 is not included in the definition of but we require it to be < 0.0155 
during the minimization, r, s, t,p and elements of H, F, G are now treated as free and the 

definition now includes the charged lepton masses as well. This results in significant 
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improvement in the fit and one is able to fit 15 observables in terms of 13 parameters with 
X = 3.01. The fit to the bottom and the down quark masses also improves. Sckm remains 
zero in this case. 



(C) For this case, we depart from the exact 23 symmetry and take /i22 different from 
/133. As already discussed, this breaks the generalized CP and results in a non-trivial CKM 
phase. Remarkably, a very small (~ 8%) breaking of the 23 symmetry is able to generate a 
non-trivial CKM phase and xLin = 3.02 with 2 degrees of freedom. Bottom quark mass is 
the only variable which deviates from its central value considerably. 



Some of the observables are not part of the their values get fixed at the minimum. 

These are shown as predictions in Table i^^. These include the CP violating Dirac phase 
SpMNS and the Majorana phases ai,2 as defined in These are trivial for the cases (A) 
and (B) due to the generalized CP invariance but one obtains non-zero values displayed in 
the Table in case (C). 



Before going into the more detailed predictions, let us underline some important points 
connected with the above fits. 



Detailed fits to fermion masses have been considered in a number of papers with 
or without 13|, ll5l] the addition of the 120-plet to the minimal 10 + 126 



19, 20 



Higgs fields. The minimal model without the 120-plet but not imposing reality of the 
coupling has more parameters than the present case but the fit is not better compared 
to here, e.g. the fit in pure type-II case [l3| with 18 parameters and 15 data points 
gives a minimum around 14.5 



The best fit solutions in cases (B) and (C) give ^23 close to maximal and close to 
zero as seen from Table (ITT 



We have fixed the overall scale of neutrino mass in eq.(l6]) by using the atmospheric 
scale as normalization. The resulting values are displayed in Table dTT])- In all three 
cases, tl comes close to 10~^°. is related to the mass of the left handed triplet resid- 
ing in the 126 representation and to other parameters in the superpotential. Detailed 
analysis 1^, ll3|, [ij, ll5|] has shown that one needs this triplet mass to be at an interme- 



diate scale ~ 10^^ GeV if the overall neutrino mass scale is to be correctly reproduced. 
The presence of such light triplet conflicts with the gauge coupling unification. An 
additional 120-plet does not qualitatively alter the situation. One possible solution 
suggested 



29 



in the literature is to add a 54-plet of Higgs and allow 150(10) to break 
first to SU{5) leaving a complete 15-plet of Higgs light at around 10^^ GeV. Other 
solution corresponds to having split supersymmetry breaking 15| . Third possibility is 
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A 


B 


C 


Quantity 


Pull 


Pull 


Pull 


md 


-1.47532 


0.167255 


0.0620115 


rris 


-0.8225 


0.271662 


-0.0545523 


nih 


-2.52388 


1.68787 


1.72811 


rriu 


0.274609 


-0.00446626 


-0.00184452 


tn 


—0 0125887 


000159604 


n 00744292 


mt 


0.00190476 


0.00901941 


-0.0199522 


rrie 





-0.000951761 


0.000179815 







0.0176266 


-0.000749102 







-0.0192274 


-0.017642 


sol 


0.679035 


-0.169337 


-0.0544521 


sin 0% 


-0.0116059 


0.00250491 


-0.00412383 


sin 


0.155231 


-0.00717926 


0.0402861 


sin 


-0.0705362 


0.0000163982 


0.0163964 


sin^ 


0.112082 


-0.111783 


-0.00578002 


sin2 0^3 





0.129873 


-0.141465 


ScKM 


— 


— 


-0.0364271 


x' 


9.80473 


3.00957 


3.02019 




Predictions 


Predictions 


Predictions 


sin^ 6*00 


0.5 






sin^ e\o 





0.000471537 


0.000226908 


^CKM 


0° 


0° 




(JpMNS 


0° 


0° 


-12.759° 


ai 


180° 


180° 


169.80° 


Q!2 


0° 


0° 


-9.445° 


rL 


2.8714 X 10"^° 


1.8183 X 10"^ 


1.8645 X 10"^ 



TABLE II: Best fit solutions for fermion masses and mixing obtained assuming the type-II seesaw 
dominance. Various observables and their pulls obtained at the minimum are shown in three cases 
(A)-(C) defined in the text. 



to allow type-I seesaw dominance 18|, |23| . We shall look at this in the next subsection 
in the present context. 

We now turn to predictions in the neutrino sector. The firm predictions of the scheme 
can be obtained by checking the variation of with the values of various observables. As 
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Parameters 


A 


B 


C 


hii 


1.95914 


-0.357916 


-0.818923 


h22 


466.637 


-649.2 


-701.354 


^23 


283.929 


-54.7552 


-32.0485 


^33 


466.637 


-649.2 


-598.783 


hi 


-1.25174 


-0.176133 


-0.343138 


hi 


14.2058 


-2.16375 


-2.07269 




-71.54 


11.5434 


11.2606 




95.5358 


-14.754 


-14.3836 




-1.66646 


3.54811 


4.19817 


523 


-26.5205 


614.356 


617.845 


r 


106.129 


61.8507 


61.1056 


s 


114.802 


-109.87 


-121.664 


t 


-1.9006 


67.0199 


65.9824 


V 


22.8456 


-0.989943 


-0.980791 



TABLE III: Values of parameters of the fermionic mass matrices in eq.(I2]) corresponding to the 
best fit solutions displayed in table([TlI). The cases(A)-(C) are defined in the text 



m 



13 



20 



we pin down a specific value po of an observable P by adding a term 



P-Po 
O.Olpo 



to then minimizing 



-2 _ 2 I 2 

X =X +Xp ■ 



If P happens to be one of the observables used in defining then its contribution is removed 
from there. Artificially introduced small error fixes the value po for P at the minimum of 
the x^. We then look at the variation of 

Xmin \X XpJlTnin 

(16) 

with pq. The results are displayed in Figs. (1-3). 



Fig.([T]) shows the variation of Xmin various pinned downed values of sin^^23- It is 
seen that the minimum occurs when sin^ is fixed to around 0.46 rather than the value 
0.5 obtained in the fits shown in Table (ITll)- The variation of Xmin is not drastic and all 
values in the range 0.3 — 0.7 are allowed at 90%CL. In comparison, variation of Xmin with 
sin^ 6'[3 shown in Fig. ([2]) is little more significant. There is a preference for values close to 
zero but values up to 0.008 cannot be ruled out at 90% confidence level. Fig. ((31) shows 
the prediction for the PMNS phase in the leptonic mixing matrix. Clear prediction is the 
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FIG. 1: Variation of Xmin with sin^ Type-II seesaw. 




0.000 



0.002 



0.004 



0.006 



sin^6i' 



0.008 



13 



FIG. 2: Variation of xLin "^i^^ 

sin^ in Type-II seesaw 



negative values for the sin 6pmns- However, all negative values are allowed within the 90% 
confidence limit. 



B. Numerical analysis: type-I seesaw 

The structure of the neutrino mass matrix in the type-I case is qualitatively different 
compared to the type-II case. Unlike Ail^ , M.^ is not /i-r invariant in general. But it can 
be made approximately /i-r symmetric if either 120 contribution or the 10 + 126 dominates 
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FIG. 3: Variation of Xmm with sytl^pmns in Type-II seesaw. 

in Mb-, see eq.([2]). We discuss below fits in three quahtatively different cases as done for 
the type-II dominance. 

(A) Here we impose the exact /i-r symmetry for by hand, i.e. by choosing g = in 
Ml). As before, \Ji is also chosen [i-t symmetric. The input parameters and observables are 
the same as in the case (A) of type-II seesaw. The results of the fits are displayed in the 
first column of the Table (lIVIl . The total involves 11 observables and is determined by 8 
parameters . The minimum value is ~ 13. While most observables can be fitted nicely, the 
top quark mass deviate by 3.6cr from the central value. Enforcing the exact /i-r symmetry 
does not appear to be a very good choice. 

(B) In this case, we do not take q = ^. A4l now satisfies eq.(l9]) and is not symmetric 
under /i-r symmetry. Ol^^ is not fixed to be maximal but is included in the definition of x^- 
As in the earlier case (B), is defined by 15 observables and is determined in terms of 14 
parameters. The CP violating phases are zero in this case and the CKM phase is therefore 
not included in x^- Experimental bound on ^^^g shown in Table ^ is imposed during the 
minimization. One now gets excellent fit to all the included variables with Xmin ~ 0.017. 

(C) In this case we introduce a small explicit /i-r symmetry breaking by assuming /122 7^ 
/i33 in eq.([2]). This allows CP violation, x^ definition now includes all 16 observables and 
depends on 15 parameters. Bound on ^'^3 is imposed during minimization. Once again we 
get an excellent fit to all the observables with Xmin — 0.18. CP violating phases in the 
PMNS matrix come as predictions. 

Noteworthy features of the fits in (B) and (C) cases above are the following: 
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A 


B 


C 


Quantity 


Pull 


Pull 


Pull 


rud 


-0.31569 


0.0346007 


-0.379829 


nis 


0.473034 


-0.0483779 


-0.0717277 


nib 


-0.108264 


-0.113763 


-0.114314 


rriu 


0.50263 


0.00026323 


0.00344698 


nric 


-0.151225 


-0.000606809 


-0.00938266 


nit 


-3.60744 


-0.0193107 


0.0122663 


nie 





-4.874 X 10"^ 


0.0000348858 







0.000480511 


0.00078371 







0.00254153 


-0.0106065 


Am? , 

sol 


-0.00977627 


-0.00192856 


0.0125218 


sin 9% 


0.0218205 


-0.00061312 


0.00761817 


sin ^no 


0.00289271 


0.00129946 


0.0284214 


sin of 2, 


-0.238953 


-0.00823361 


0.0366413 


sin^ 6*^2 


-0.0129712 


0.000590904 


-0.00265193 


sin2 





-0.00544523 


0.0289959 


^CKM 


— 


— 


-0.120278 




13.6821 


0.0169632 


0.180526 




Predictions 


Predictions 


Predictions 


sin^ 


0.5 






sin^ 6lio 





0.0135605 


0.013505 




0° 


0° 




^PMNS 


0° 


0° 


-0.287748° 


ai 


180° 


0° 


2.156° 




0° 


0° 


2.616° 


rR 


4.1143 X 10"^^ 


5.2329 X 10"^^ 


5.0093 X 10"^^ 



TABLE IV: Best fit solutions for fermion masses and mixing obtained assuming the type-I seesaw 
dominance. Various observables and their pull obtained at the minimum are shown in three cases 
(A)-(C) defined in the text. 



The overall neutrino mass scale is determined to be around ~ 5 x 10^^^. rpi 
is related to the ratio of the vev of the doublet and the RH triplet components 
in 126. The values of r^j obtained here are similar to the values obtained in 
201] which assume 126 RH triplet vev to be at the GUT scale. Thus one does 
not need an intermediate scale in order to fit the neutrino masses and one can ob- 
tain the gauge coupling unification. This is consistent with observations in 
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20 



23 



14 



• Maximality of ^23 is not imposed. But it is fixed to be very close to j at the minimum 
in both the cases. The departure from the /i-r symmetry results in being non-zero 
and is fixed around the upper bound at the minimum as seen from the Table flIT|) . 



Although an explicit breaking of the /x-r symmetry is introduced in case (C), the 
amount of the breaking required in order to obtain the large CP violating phase is 
extremely tiny, 

0.0045 . (17) 



^22 ~ ^33 



^22 + ^33 



The exact fi-r symmetry is known [5[ to lead to the unwanted predictions Vub = Vcb = 
sin^ 6*23 = 0. Here we have two sources of breaking this symmetry, spontaneous through 
the vev of the 120-plet and explicit through eq. (fT7|) which allows one to reproduce the 
mixing angles correctly. In spite of the /x-r breaking, the final fermion mass matrices 
display a remarkably good fi-r symmetry. We make this explicit by giving the quark 
and lepton mass matrices in the case (C) above in Appendix (A). Mu^d,i and Ail are 
seen to be nearly fi-r symmetric. There is an order of magnitude difference in the 
imaginary parts of the 12 and 13 elements of Adl- But these imaginary parts are 
much smaller than the corresponding fi-r symmetric real parts. The only source of 
the large /i-r breaking occurs as a difference between the 12 and 13 elements of the 
Dirac neutrino mass matrix Mo. This results from the spontaneous breakdown and 
rather large value of the parameter q. 



As in |19l . |20| we have concentrated here in obtaining generic fits to fermion masses 
rather than considering the entire parameter space of the theory given by the Yukawa 
couplings and basic parameters in the superpotential. Parameters in fermion mass 
matrices are related to the strengths of the light Higgs components in various S'C'(IO) 
Higgs representations. These are determined by the fine tuning conditions and the full 



superpotential. Grimus and Kiihbock [l9j have laid down consistency constraints on 



these parameters following from these fine tuning relations and from the requirement 
that the Yukawa couplings stay in the perturbative regime. We have checked that 
these conditions are satisfied by the parameters given in the Table fllllfVp . 



We follow a similar procedure as in the type-II case to obtain possible predictions on the 
neutrino mixing variables. We pin down an observable P to a specific value po by adding 
a contribution Xp to x^- We then determine the variation of Xmin defined earher with 
Po- Variations of Xmin obtained at different local minima are shown as scattered plots in 
Fig. (4-6). 



Clear predictions emerge unlike in the type-II case. As Fig.Q shows, the sin^ ft 



23 



IS 



preferentially restricted near 0.5 and one obtains the limit ~ 0.42 — 0.63 at the 90% CL 

^3 



Fig.© shows similar variation with respect to sin^^^3. Here, the preferred values occur 
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Parameters 


A 


B 


C 


/ill 


907.294 


34.9749 


35.0178 


/l22 


119.541 


554.305 


556.777 


^23 


-119.052 


554.457 


554.429 


^33 


119.541 


554.305 


551.775 


/ll 


74.5214 


-15.7284 


-15.716 


/l2 


-2.82327 


20.8852 


20.8951 


f22 


-74.237 


-29.4577 


-29.4636 


/23 


74.2104 


-29.5265 


-29.5305 


5l2 


182.676 


3.10944 


2.79728 


523 


-4.5309 


-3.5854 


-3.21385 


r 


1.24579 


83.0642 


83.7973 




0.266298 


176.883 


178.571 




0.844656 


0.450978 


1.0715 




2.35413 


0.0117737 


0.011244 


Q 





4042.93 


4537.34 



TABLE V: Values of parameters of the fermionic mass matrices in eq.([2]) corresponding to the best 
fit solutions displayed in table (llVjl . The cases (A)-(C) are defined in the text 



near the present limit and one obtains sin > 0.005 at 90% CL. The predicted values for 
smSpMNS are displayed in Fig. ([6]). These are negative but very small. 

All the above solutions are obtained through an extensive search using the random search 
algorithm in Mathematica and the MINUIT sub-routine in FORTRAN and we have shown 
in tables the solutions corresponding to the minimum that we obtained. Considering the 
non-linearity and complexity of the problem here, it is difficult to rule out the existence of 
still lower minima and predictions may improve if they exist. 

We end this section with a comment on the specific /i-r symmetry defined by S used in 
eq.(l5]). Definition of S is basis-dependent. One could change the original basis of the 16-plet 
through an arbitrary rotation R. The structure of the Yukawa couplings and the resulting 
fermionic mass matrices would look different in the new basis. The new Yukawa couplings 
would still satisfy the same equation as (jlj) but now with a rotated 5*: Sr = R^SR. Thus 
the /i-r symmetry may appear to look different with different choices of R. Specifically, if 
R corresponds to a rotation by ^ in the 23 plane then the Sr assumes the form 

/l o\ 



-1 
1 



(18) 



This is nothing but the Z2 symmetry imposed in 19|] which is thus equivalent to the gener- 
alized /u-r symmetry considered here if both remain unbroken. Difference arises after these 
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FIG. 4: Variation of xtun with sin^ Bin in Type-I seesaw. 
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FIG. 5: Variation of Xmin with sin^ 6[r^ in Type-I seesaw. 



symmetries are broken. Ref. ISj] uses complex vev to achieve Z2 breaking as a result of 
which analogue of eqs.([71l8]) do not hold in their case. In our approach, we introduce small 
explicit breaking of /i-r symmetry in H. The model in 19|] has 20 free parameters compared 
to 15 used here. 

Note that the explicit breaking of the /i-r symmetry is technically natural in the super- 
symmetric context. Alternatively, one can achieve such breaking by introducing an addi- 
tional 10-plet of the Higgs field which changes sign under the /i-r symmetry. Combined 
contributions of these two 10-plets would then give an exphcitly /i-r non-invariant H. 
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FIG. 6: Variation of Xmin with sin(5pAfArs in Type-I seesaw. 



IV. SUMMARY 

Aim of this paper was to integrate the successful /x-r symmetry within the SO{10) frame- 
work in order to obtain a constrained picture of fermion masses and theoretical understand- 
ing of the largeness of the atmospheric mixing angle. The explicit model discussed here pro- 
vides this integration rather well as shown by the detailed fits to fermion masses presented 
in Tables (]II|IV|) . Interestingly, mass matrices obtained in the model under consideration 
display a generalized CP invariance if Yukawa couplings are taken to be /i-r symmetric. 
Small explicit breaking of this symmetry is sufficient to generate the required CP violating 
phase. The best scenario is obtained in the type-I seesaw model with very tiny explicit fi-r 
symmetry breaking. This scenario is characterized by the predictions sin^ ~ 0-42 — 0.63, 
sin^ > 0.005 and negligible CP violation in neutrino oscillations. Final quark, the charged 
lepton and the light neutrino mass matrices (collected in the Appendix) respect /x-r symme- 
try to a very good approximation indicating that this symmetry provides a good description 
of the entire fermion spectrum rather than being restricted to the neutrino sector alone. 
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VI. APPENDIX (A) 

We list here the fermion mass matrices following from eq.([2]) using the best fit values 
of the parameters given in Table (jV]) corresponding to the type-I seesaw mechanism. The 
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neutrino mass matrix is expressed in eV units while all other mass matrices are expressed 
in MeV units. 



/ 19.3018 20.8951 + 2.79728i 20.8951 - 2.79728i 

20.8951 - 2.79728i 527.314 524.898 - 3.21385i 

\ 20.8951 + 2.79728i 524.898 + 3.21385i 522.311 

/ 127.971 3731.25 + 2. 99727i 3731.25 - 2. 99727i 

3731.25 - 2.99727i 41395.1 41186.3 - 3.44363i 

V 3731.25 + 2.99727i 41186.3 + 3.44363i 40975.9 



M, 



Md 



82.1659 -62.6852 + 0.0314526i -62.6852 - 0.0314526i 

-62.6852 - 0.0314526i 645.168 643.02 - 0.0361365i 

-62.6852 + 0.0314526i 643.02 + 0.0361365i 640.166 

/ 11353.7 -11193.7 + 12692. 2i -11193.7 - 12692. 2i 

= -11193.7 - 12692. 2i 62440.4 62279.4 - 14582. 3i 

\ -11193.7 + 12692. 2i 62279.4 + 14582.3i 62021.3 



(19) 



(20) 



(21) 



(22) 



Mi 



-0.0242264 
-0.0143681 + 0.0004742i 
-0.0143657 - 0.0000755678^ 



-0.0143681 + 0.0004742i 
-0.0128288 + 0.00678282i 
-0.0163109 + 0.000214216i 



-0.0143657 - 0.0000755678i 

-0.0163109 + 0.000214216i 

-0.0127693 - 0.00629525i 

(23) 
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